Disorder in quantum vacuum: Casimir-induced localization of matter waves 
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Disordered geometrical boundaries such as rough surfaces induce important modifications to the 
mode spectrum of the electromagnetic quantum vacuum. In analogy to Anderson localization of 
waves induced by a random potential, here we show that the Casimir- Polder interaction between a 
cold atomic sample and a rough surface also produces localization phenomena. These effects, that 
represent a macroscopic manifestation of disorder in quantum vacuum, should be observable with 
Bose-Einstein condensates expanding in proximity of rough surfaces. 
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Introduction- Waves propagating in disordered poten- 
tials undergo multiple scattering processes that strongly 
affect their usual diffusive transport and can result in 
localized states. In one dimensional systems Ander- 
son localization is a ubiquitous phenomenon jl^. Re- 
cently it has been observed in a ID Bose-Einstein con- 
densate (BEC) expanding in the presence of a random 
laser speckle field [2| or a bi-chromatic potential [3] . The 
asymptotic density profile shows exponential localization 
even in the weak disorder limit [^]. In this Letter we 
show that the high sensitivity to disorder of ID cold 
atomic systems can be strong enough to yield localiza- 
tion of matter- waves due to disorder in vacuum. 

The Casimir-Polder (CP) interaction between a BEC 
and a flat surface has been recently measured (see [5] and 
references therein). Cold atoms act as local probes of the 
electromagnetic quantum vacuum, and they have been 
proposed as good probes of the influence of non-trivial 
geometrical effects on the CP atom-surface interaction 
potential [6|, [D]. Charge and current quantum fluctua- 
tions in a rough surface induce a disordered CP poten- 
tial that directly affects the dynamics of a BEC trapped 
close to such a surface, and may lead to localization of 
matter waves. In this Letter we will investigate the lo- 
calization properties of the BEC, which reveal the CP 
disordered interaction with a random surface. An es- 
sential tool for computing the Casimir-Polder potential 
above a disordered surface is the scattering approach [8| , 
that computes the Casimir interaction between bodies as 
a non-trivial multi-scattering process. 

Atom- surface Casimir interaction- For simplicity, we 
will consider a surface with translational invariance along 
some direction (say ?/), and a generic uni-axial profile 
h{x) = 5^^-,^ /li cos(/ciX + ^i) (see Fig. 1). Here hi are 
the amplitudes of the Fourier spectrum of the profile, 
\i = 27r/ki the corresponding periods, and 6i are off- 
sets. In the case of a corrugated surface these are fixed 



parameters, while for stochastic roughness they are ran- 
dom parameters distributed according to certain proba- 
bility distributions. In the following we will consider the 
latter case, with flat and independent probability distri- 
butions in certain intervals (which is the simplest noise 
model to describe stochastic surfaces). A ground-state 
atom above such a surface is affected by a CP potential 
/7(x,z), where z is the atom-surface distance (defined 
from a mean surface). The exact computation of the 
CP potential for certain corrugated profiles, such as ID 
lamellar gratings [9|, llO| , can be performed via scattering 
theory. For stochastic roughness, however, exact results 
do not exist. Previous works have computed the CP po- 
tential U{x, z) by means of the pair-wise approximation 
ll| . Here we calculate it using the scattering approach, 
that takes into account non-additivity effects. We ex- 
pand the potential to second order in powers of /i^, that 
we assume are the smallest length scales in the problem. 
In general, U{x, z) has a part that depends only on z 
(that for planar surfaces gives the usual CP force) and 
one that depends both on x and z, that gives a lateral 
component to the CP force due to the lack of transla- 
tional invariance. We shall denote this latter component 
as Ul{x^ z). The first order was already computed in [6[, 

U^l\x,z) = -^^f2h,g^'\k,z) cos{k,x+0i), (1) 

where g'^^^ is a response function that depends on the op- 
tical properties of the surface. Here we need the second- 
order correction [12 1 
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FIG. 1: a) Cigar-shaped BEC of width a and axial size 2Ltf 
above rough surface. Here Ltf is the initial Thomas- Fermi 
length of the BEC. Dimensionless response functions for a 
perfectly reflecting surface: b) g^^'{Z^ (^ Z — kz) and c) 
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^'^\Zi^Zj) ( Zi — kiZ and Zj = kjZ.) 



where the response function ^^^^ also depends on the 
material properties. In Fig. 1 we show g^^"^ and ^^^^ 
for a perfectly reflecting surface. The dimensionless ker- 
nel g^^\kz) decays exponentially for large kz (e.g., for 
a perfectly reflecting surface it is g^^\Z) = e~^(l + 
Z + 16.EV45 + ^^/45), with Z = kz). The dimension- 
less kernel g^'^\kiZ^kjz) typically decreases for large z 
along generic directions in the ki^kj plane, but it grows 
along the ki = —kj direction. However, since there is 
a maximum k admissible in the perturbative expansion 
(roughly given by 27r//i, where h is the typical order of 
magnitude of the height profile), this imposes a maxi- 
mum value to ki^kj given by 27r//i, and therefore there 
exists an upper bound to the response function given by 
g^'^\27rz/h^—27Tz/h). Taking this fact into account, it 
follows that the second order perturbative expansion is 
sufficient for the purpose of this work [12]. 

Localization due to Casimir- Polder interactions- We 
consider the expansion of a tightly confined cigar-shaped 
BEC parallel to the rough uni-axial surface, so that 
its axis is perpendicular to the corrugation lines (Fig. 
1). The effective ID dynamics of the dilute BEC can 
be described by a mean- field wavefunction (p{x^t) which 
evolves according to the ID Gross-Pitaevskii (GP) equa- 
tion [135 : 

ihdt(p{x,t) = --— S^(/:)(x,t) + UL{x,zo)0{t)(p{x,t) 
2m 



muj'ix'^ 



Q{-t) (p{x, t) + ge^\^{x, t)\^if{x, t), (3) 



where ^eff = ^/27rcr^ is the effective coupling constant 
for the ID problem (a is the width of the radial Gaus- 
sian profile). Although CP interactions are known to 
be non-additive, for dilute BECs additivity is a good 
approximation [1^. We assume that the harmonic po- 



tential muj^x^ /2 confines the system for t < 0, so that 
(f{x^t) = e~*^*/^(po(^) with /i the chemical potential. 
At t = the axial trap is turned off causing the BEC 
to expand in the presence of the disordered potential 
Ul{x^zo). Note that the distance to the surface, zq, in 
Eq. (|3|) is fixed during the evolution (the tight confine- 
ment freezes the radial motion), that is why all constant 
potential terms, including the x- independent component 
of the CP potential, are absorbed in a global phase-shift 
of the GP wavefunction. We intend to study the im- 
print left by CP forces on the asymptotic density pro- 
file n{x^t) = \(p{x^t)\'^ of the BEC after the expansion 
through the disordered potential. 

More specifically we study the BEC density profile av- 
eraged over many realizations of h{x)^ that we denote as 
n{x). One of the parameters characterizing the evolution 
of the BEC is the root mean square Vr of the disordered 

potential, defined as V-^{zo) = {Ul{x, zq) — Ul{x, zq))'^ 
{x dependence disappears due to the random offsets Oi). 
When the strength of the disorder Vr{zo) is much smaller 
than the chemical potential, the evolution can be treated 
analytically in perturbation theory [^ . The applicability 
of such perturbative approach is vahd in the CP context 
when relatively large distances zq are considered (note 
that lim^Q^oo "^(^o) = 0)- However, when Vr{zo) ~ /i 
(small values of 2^0), a full numerical approach is required. 
Let us first consider the weak disorder case, Vr(^o) ^ 
II. When the trap is switched off the first stage of the ex- 
pansion is driven by interaction term gQ^\if{x.,t)\^ip{x.,t) 
in (|3]) and the disordered potential can be neglected. 
After the initial expansion process the density becomes 
low enough to make the disordered term dominant and 
thus the system can be considered as a collection of 
non-interacting particles expanding in a weak disor- 
dered potential. In this situation the averaged den- 
sity n{x) of localized atoms is fully characterized by 
the two-point correlator of t he disordered CP p oten- 
tial [4]: C(\x - x'\]Zq) = U l{x , zq) U l{x' , zq) . In 
the following we will omit writing the explicit depen- 
dence on the parameter zq. This quantity determines 
to which extent a non-interacting particle with momen- 
tum k localizes in the random potential. The informa- 
tion about the localization length is given by the in- 
verse of the so-called Lyapunov exponent jl5|, 7(/c) = 
{m/Ah^Ek) /^ C{x) cos{2kx) dx, where Ek = h^P/2m. 
The consistency of this perturbative approach requires 
7(/c) <C /c to hold, that is the typical length scale of the 
exponential localization must be larger than the typi- 
cal wave-length of the unperturbed single-particle state. 
Deviations from a pure exponential behavior were calcu- 
lated in 16| , and used in [J] to find an expression for the 
averaged density profile, namely: 

ee)\M^Wdk, 
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where ^ = h/ ^Amfi is the healing length of the BEC 



and the averaged modulus squared of the single particle 
wavefunction (j)^ is 
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Using these relations n{x) can be calculated once C{x) 
is given. In general, two different asymptotic regimes for 
the wings of the density profile can be identified, depend- 
ing on the healing length of the BEC, the spectral content 
of C, and the maximum length scale Lmax that can be 
measured in the system (e.g. the spatial observation win- 
dow of the measurement set up). When the above system 
is solved for n{x) assuming that 7(/c) 7^ OV/c G (0,1/^) 
and Lmax ^ 00 an exponentially localized profile is al- 
ways found. On the other hand, when ^(k) vanishes in 
an interval k < k < 1/^ and Lmax -^ 00, there exists 
an effective sharp mobility edge at /c = /c, the modes 
with k > k continue to expand and the localized modes 
with k < k give rise to an algebraic-decay of the form 
n{x) oc 1/x^. This is typically the case of the correlator 
of the laser speckle potential, that has a sharp cut-off. In 
contrast, the function j{k) for the CP disordered poten- 
tial of the stochastic surface considered here reads: 



j{k) = 



.277^2 
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h{x)^P{7r/k) (^g^'\2kzo))\ (6) 



where h{x)'^ is the mean of h{x)'^ at any point x, and P(A) 
is the probability density for each A^ (in this case flat) 
and F{zo) = 3fica(0)/87r^eo^o- Hence 7 is exponentially 
suppressed (see Fig. [2]). In the CP case exponential local- 
ization is also found if Lmax -^ 00. However, for a finite 
value of Lmax7 the situation may change. In fact, in typi- 
cal cases where the perturbative approach is valid in this 
context the Lyapunov exponent satisfies 7(l/0^max <^ 1 



making the modes with k^ 



> 



1 



"Hl/^max) effectively de- 



localized, see the inset in Fig. [2] (here 7"^ denotes the 
inverse function of 7). Thus, in a restricted region of size 
2i^max, the function n{x) does not decay exponentially 
but rather algebraically (this assertion is approximate 
and becomes exact in the limit 7(l/0^max -^ 0), the 
fast decay of j{k) results in an effective mobility edge 
at /c*. However, the transition between the two limiting 
regimes is not sharp in the CP context, and in principle 
one can pass continuously between both varying Lmax- 
Finally, note that this is an external parameter which 
does not affect the system but only has to do with the 
measurement process. 

We now compare this perturbative approach with a full 
numerical simulation. We consider a BEC of iV = 10^ 
^^Rb atoms tightly confined in the radial direction. The 
width of the ground state is assumed to be a = 0.25/im 
(i.e., radial trapping frequency o;^ = 27r x 286Hz) and 
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FIG. 2: 7(/c) vs. k^ from Eq.(|6|. The atom-surface distance 
is zo = 1.5/im. Inset: same data in Log-Lin scale, the maxi- 
mum length scale to be measured Lmax = 1mm, and the value 
/c* = 7~^(l/i^max) separating localized (Loc.) from delocal- 
ized (Deloc.) modes. 




FIG. 3: BEC density (arbitrary units) vs. position. Both 
the perturbative theory described by Eas. (|4l5l6[) (solid) and 
the full numerical simulation (dots) are computed using the 
first order approximation for the CP potential at zo = 1.5//m. 
The surface profile is averaged over 40 realizations. Time 
corresponds to cj^t = 28. 



its axial size 35/im (i.e., c^a^ = 27r x 2.75Hz). For these 
values the healing length of the BEC is ^ = 0.85/im. A 
perfectly conducting uni-axial rough surface is brought 
to close proximity of the BEC, with typical distances of 
Zo ^ 1/im. The stochastic surfaces are generated using 
between 15 and 25 harmonics, and the parameters /i^, A^, 
and Oi are taken as independent random variables with 
flat probability distributions satisfying hi G [0,200]nm, 
6i G [0,27r], A^ G [Amin,Amax]- Notc that the approxi- 
mation we are using to evaluate the CP lateral poten- 
tial assumes Amin > h^ however this is not an impor- 
tant restriction because the exponential suppression of 
the modes with the factor 2^0 /Amin makes the evolution 
insensitive to the lower limit. In Fig. 3 we compare a di- 
rect numerical simulation with the perturbative approach 
for the function j{k) shown in Fig. 2, which corresponds 
to a surface with 25 harmonics and A^ G [l,20]/im at a 
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FIG. 4: BEG density (arbitrary units) vs. position, after 
Uxt = 14. The profile is averaged over 40 realizations at a 
distance of zq = 1.0/im. Numerical simulation (dots) includes 
both the first and second order terms of the later al GP poten- 
tial Ul{x, z). The wing is fitted by a power law n{x) oc 1/x^ 
with v — 1.84 (solid). Inset: zoom of the numerical data in 
Log-Lin scale. 



distance of zq = 1.5/im of the BEG (in this situation the 
mean square root of the noisy potential is l^ = 0.089 ji). 
Even in this pertubative regime the integral in ([5]) must 
be calculated numerically. In the preceding case we have 
taken into account only the first order correction to the 
lateral CP potential because it is a good approximation 
in such regime (this was verified in the numerical simu- 
lations). Note t hat, wh ile the BEG evolves, the averaged 
density profile n(x,t) (here truncated after 40 realiza- 
tions of the surface profile) approaches to the asymptotic 
t ^ oo solid line n{x) predicted by the theory [4]. Be- 
cause of the finite value of the time variable used here, 
uoxt = 28, it can be seen that the edge at the right is 
expanding yet, snapshots at shorter times show a bigger 
elbow for smaller positions. The space window chosen 
corresponds to L^ax = 1mm and, in that scale, a log-log 
graph makes the prediction of the perturbative approach 
to be a straight hue (slope constant within 1% accuracy). 

In the following we consider stronger perturbations, so 
that Vr{zo) is no longer much smaller than /i. In such cir- 
cumstances the previous perturbative treatment breaks 
down since the GP potential cannot be disregarded in the 
first stage of the expansion. Therefore a fully numerical 
method is required to solve Eq. (3). For distances of 
about zo ^ 1/im the second order correction to the lat- 
eral GP potential Eq. (J2j) also becomes important, so in 
the full numerical simulation we solve exactly the evolu- 
tion of the BEG taking into account these second order 
effects. The results are shown in Fig. 4, where we plot a 
typical stationary density profile averaged over 40 real- 
izations of the surface at zq = 1/im from a surface with 15 
harmonics in the range A^ G [1, 8]/im and hi G [0, 200] nm 
(for these parameters Vr = 0.5/i). Most of the atoms 
are trapped in a high density core of typical size given 



by the initial size of the BEG due to single barrier re- 
flections (similar to the strong disordered regime consid- 
ered in y ) , the atoms leaving the core develop a density 
profile with algebraic-like wings. Again the exponential 
suppression of the GP potential makes the evolution in- 
sensitive to modifications in the lower limit of the range 
of wavelengths Amin, the spatial cut off frequency be- 
ing imposed by zq (this was also confirmed numerically 
pushing the surface spectrum to shorter wavelengths). 
There is a high sensitivity with the distance to the sur- 
face due to two effects, the strong power-law decay in zq 
of the GP potential, and the exponential suppression of 
the modes with kzQ. Their combined effect makes the 
system pass rapidly from the perturbative to the non- 
perturbative regimes. In fact it can also be seen that for 
configurations such as the one presented here distances 
of about Zq ~ 3/im make Vr so small compared with the 
typical kinetic energy that the effect of the GP potential 
on the expansion of the BEG is negligible. Finally, it is 
worth mentioning that the localization phenomena found 
for the parameters used in our work may break down for 
stronger BEG non-linearities ^eff in Eq. (jSj jlTj . 

We have shown how quantum fluctuating interactions 
in non-trivial geometries affect transport properties of 
matter waves. We have found that a quasi- ID Bose- 
Einstein condensate expanding close to a rough surface 
undergoes algebraic localization. The perturbative pre- 
diction for the density profile of the BEG localized via 
the disordered Gasimir-Polder interaction was shown to 
be in good agreement with a fully numerical approach. 
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